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Abstract 

A search for J ftp radiative decay to weakly interacting neutral final states was performed using 
the CLEO-c detector at CESR. J/ip events were selected by observing the hadronic decay ip(2S) — > 
7r + 7r~ J/ip. A total of 3.7 x 10 6 J/ip events was used to study the decay J/ip — > 7 + X, where X 
is a narrow state that is invisible to the detector. No significant signal was observed and upper 
limits on the branching fraction were set for masses mx up to 960 MeV/c 2 . The upper limit 
corresponding to mx = is 4.3 x 1CT 6 at the 90% confidence level. 
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Understanding the nature of dark matter is one of the major goals of particle physics. In 
one scenario two dark-matter particles annihilate in the early universe primarily through the 
s-channel production of a light neutral boson that lies outside the standard model. Coupling 
strengths can be chosen that allow both the dark-matter particle and the boson to decouple 
from standard-model particles and thus avoid conflict with experimental limits on their 
laboratory production rates, while at the same time giving the correct remnant dark- matter 
density [lH8|. In these models the lightest neutralino is a logical candidate for the dark- 
matter particle but a light scalar is also a possibility. Non-minimal supersymmetric models 
that include a scalar singlet Higgs field produce a rich spectrum of Higgs and neutralino 
particles that can extend to very low masses 0, 1(| so they can be easily applied to the 
dark-matter problem. For example, the required light boson can be a gauge boson 04^, @] 
or a CP-odd Higgs boson [l|, @] . 

Experimental limits on the branching fractions for quarkonium radiative decay to invisible 
final states (H, 12] provide laboratory constraints for some of the models discussed above 0, 
0]. Preliminary data on T(35) — > 7 + invisible have also been reported [l3j]. These can be 
incorporated into theoretical models by calculating the branching fraction for a light boson 
to be produced in radiative quarkonium decay and then assuming its dominant decay is to 
two neutralinos. The neutralinos interact very weakly with normal matter so they would 
not be detected in an experiment. In this case the charmonium radiative decay amplitude 
takes the form, 

T(J/VW7X) _ G F m 2 c cos 2 6 



r( J /if; ->■ fi+fi-) 2V2air tan 2 (3 

where X is a boson with pseudoscalar quark coupling, 9 is the Higgs mixing angle, m c is half 
the J /if) mass, and tan/3 is the usual ratio of vacuum expectation values. A factor of 1/2 is 
included in this expression as a rough estimate of radiative and relativistic corrections j3]. 
The corresponding expression for bottomonium decay takes the same form but with tan 2 
appearing in the numerator instead of the denominator (Eq. (Ill) in Ref. [H). 

In this letter we report upper limits on the branching fraction for J/tp — > 7 + invisible. 
Here an invisible particle is one that escapes detection because its interaction with the 
detector is very weak, or because it decays to other neutral weakly-interacting particles 
within the detector volume. The present data provide new constraints for some models in 
which a light boson decays to two dark-matter particles. 

Data were acquired at thejb(2S) mass using the Cornell Electron Storage Ring (CESR) 
with the CLEO-c detector [lj-ll6|. Photons were detected in a CsI(Tl) electromagnetic 
calorimeter, which has photon energy resolution equal to 2.2% at 1 GeV and 5% at 100 MeV. 
For measurements in the charm region the CLEO III silicon vertex detector [l7[ was replaced 



with a cylindrical drift-chamber and the solenoid magnetic field was set to 1.0 T. 

J/tp events were tagged by measuring the charged pions from ip(2S) — > tt + tt~J/iIj. All 
events that had additional charged tracks beyond the tagging pions were rejected. Event 
selection criteria were chosen to optimize the identification of charged pions and minimize 
the background from hadronic showers produced in the calorimeter. Transition pions were 
required to have differential energy loss signatures consistent with a pion and their trajec- 
tories were required to form a single decay vertex close to the e + e~ interaction point. In 
addition, individual pion momentum vectors were required to lie in the central region of 
the detector with |cos#'| < 0.83, where 9' is measured relative to the beam axis, and the 
summed momentum vector for each pion pair was required to have |cos#'| < 0.95. The lat- 
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FIG. 1: Recoil mass calculated from the di-pion four- momentum, after pion selection. A fit to the 
unshaded region of the spectrum was used to extract the tagged J/ip yield. The solid line shows 
the results of the fit and the dotted line shows the background contribution. Data in the shaded 
regions were used to study background. 

ter constraint was used to suppress background from 77 fusion and direct J/tp production 
from e + e~ — > jJ/ip- Figure [1] shows the recoiling invariant mass calculated from the four- 
momentum of the tagging pions and the center-of-mass energy of the initial e + e~ system. 
This recoil mass was required to be within ±5 MeV/c 2 of the J/ip mass for the selected 
invisible-decay candidates. Further background reduction was achieved by requiring the in- 
variant mass of the 7r + 7r~ pair to be between 460 and 590 MeV/c 2 , and requiring charged 
tracks to have momentum component transverse to the beam in excess of 100 MeV/c. Those 
criteria were determined by maximizing S 2 /B, where S is the inclusive signal yield in the 
J /if) mass peak and B is the scaled background yield in dipion recoil mass side-band regions 
(Fig. [p. The tagged J /ip yield, 3.7 x 10 6 , was determined by fitting the data in Fig. [T]to a 
sum of three Gaussians and a third-order polynomial background function. 

The most energetic calorimeter shower in each event that was not associated with a 
transition-pion track was designated as the signal photon corresponding to J/ip radiative 
decay. This signal photon must lie in the barrel region of the detector, with |cos#'| < 0.79, 
and when transformed into the J/ip rest frame it must have energy exceeding 1.25 GeV. The 
latter selection avoided the rapidly rising background below 1.25 GeV. The photon shower 
was also required to have a lateral shape consistent with that expected for a photon so as 
to suppress anti-neutron showers from the decay J/ijj — > nn. 

Because the signal-event topology includes only transition pions and a signal photon vis- 
ible in the detector, restricting additional calorimeter activity should in principle only aid in 
rejecting background from other J/ip decays and not affect the invisible signal. However the 
transition pions themselves interact in the calorimeter and produce shower fragments. Hence 
care was taken so as not to reject events in which such remnants of hadronic interactions 
are produced, while still eliminating background events with additional photons. Showers 
were considered as photon candidates only if they were not in close geometrical proximity 
to the calorimeter entry point of a transition pion, and if they survived both a lateral-profile 
selection and the application of a neural network. The lateral-profile restriction eliminated 
broad showers from consideration, and the neural network examined the energy distribution 
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among the crystals in a shower. (See Fig. 12 in Ref. [18J.) Events were removed from 
consideration if there were any additional photon candidates with energy greater than 50 
MeV. This restriction was found to reject just 1% of a test sample selected by requiring 
J '/if) — >■ decay instead of the signal-photon requirement. 

In addition to the photon selection criteria discussed above, a final selection was made 
on the total energy deposited in the calorimeter minus the signal-photon energy. Events 
above 600 MeV were rejected, which reduced the background from events in which one or 
more energetic photons coincidentally overlapped the shower caused by one of the transition 
pions. This removed much of the high-energy background in the signal-photon spectrum 
arising from radiative J/if> decay to n and rj . The limit was set by comparing the results 
of Monte Carlo simulations for these two channels and for radiative decay to a massless 
invisible particle. The signal-photon selection efficiency, e, fell from about 65% to 55% after 
this data selection was made. 

Figure EH a ) shows a histogram of the selected data for J /if) — > 7 + invisible as a function 
of E*, the photon energy in the J /if} rest frame. A total of 73 data events were measured 
with E* > 1.25 GeV. A sharp peak in this spectrum would be evidence for J /if) radiative 
decay to a narrow final state. No obvious peak corresponding to J /if) — > 7 + invisible is 
apparent in the data. A Monte Carlo simulation of the background spectrum, normalized 
to integrated beam luminosity, is also shown. 

To study the background from channels such as if)(2S) decay, data were extracted from 
recoil-mass regions adjacent to the J /if) (Fig. [p. In Fig. [2] this sideband yield has been 
scaled down by the ratio of the mass intervals used to select the on-resonance and sideband 
data. The on-resonance data are seen to originate predominantly from J /if) decays, as the 
sideband yield is negligibly small. 

Although there are no obvious peaks in the on-resonance data, there does appear to be a 
relatively smooth background, the size of which is underestimated by the Monte Carlo sim- 
ulation. This background has the somewhat surprising property that it continues well above 
the photon kinematic limit of E* = 1.55 GeV, meaning it cannot be solely from radiative 
J /if) decays. Monte Carlo studies indicate that the dominant background contribution in 
Fig. |2J a), in fact the only one for E* > 1.4 GeV, is from J/if> — > nn decay. These events 
have the neutron leaving no detectable signature in the tracking or shower detectors, and 
the anti-neutron producing a shower in the calorimeter that is misidentified as a high-energy 
photon. Since anti-neutrons annihilate in the calorimeter they can produce shower energies 
that exceed the kinematic limit for photons by converting part of a nucleon rest mass into 
visible energy. This background source was studied by selecting calorimeter showers with 
large lateral spread, effectively eliminating real photon showers and preferentially selecting 
background events that did not have a photon in the final state. The E* distribution of 
these events appears in Fig. Mjo) for data and Monte Carlo, where two salient features can 
be observed: it shows no cutoff at the kinematic limit for photons, just as in Fig. EH a), an d 
it is a smooth function of E* The GEANT3 Monte Carlo simulation [l9[ correctly models 
the spectral shape but not the absolute number of anti-neutrons that satisfy the photon-like 
lateral-shower criterion. However since this is a very small fraction of all anti-neutrons, it is 
not unexpected. 

In order to extract upper limits on the branching fraction for J /if) — > 7 + invisible a 
series of fits was made to the selected on-resonance data. Two-body radiative J / if) decay to 
a narrow final state features a monochromatic photon in the J /if) center-of-mass frame. This 
produces a peak in the measured E* spectrum with a line shape determined by the calorime- 
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FIG. 2: (a) Invisible-decay events as a function of photon energy in the J/ip rest frame, and (b) 
the anti-neutron-enhanced spectrum obtained by selecting calorimeter showers with large lateral 
extent. The solid points are data and the unshaded histograms show the results of Monte Carlo 
simulations. The shaded histograms show the small contributions from events that have di-pion 
recoil mass adjacent to the J /ip mass window. 




ter resolution. We extracted the branching fraction for J/ip — >■ 7X for fixed mass, mx, with 
a binned maximum likelihood fit to the E* spectrum. The fit assumed an exponentially 
falling background and a peak lineshape determined from a Monte Carlo simulation of the 
detector response to signal events. The peak lineshape was parameterized with a Crystal 



Ball function [20j, |21|, which features a low-side tail joined smoothly to a Gaussian core. 
The E* peak resolution was approximately 30 MeV. 

Fits were performed on the data in the E* range from 1.25 to 1.65 GeV, with fixed peak 
energies chosen from 1.400 GeV to 1.548 GeV in 5 MeV steps. A typical fit is shown in 
Fig. [31 This region encompasses the kinematic threshold at 1.55 GeV and probes a broad 
range of light invisible-particle masses recoiling against the photon, up to 960 MeV/c 2 . The 
data in this range are well described by an exponential background curve alone, showing 
no evidence for a signal. We also note the possible presence of a narrow structure in the 
spectrum at about 1.36 GeV. This peak is much narrower than the expected photon-energy 
resolution, and is therefore attributed to a statistical fluctuation in the data. 

The branching fraction for invisible decay was calculated as B(J/ip — > 7 + invisible) = 
N inv /(eN tag ). Here N tag is the number of tagged inclusive J/ip events, N inv is the number of 
fitted events after photon selection and particle vetoes, and e is the efficiency for selecting 
invisible radiative decays from the tagged J/ip sample, as determined from Monte Carlo 
simulations. 

As there were no strong peaks in any of the fits, branching fraction upper limits were 
extracted at the 90% confidence level at each mx- This was done by integrating the likeli- 
hood for positive branching fractions. The limits are dominated by statistical uncertainties 
but vary somewhat depending on the shape of the background function (exponential or 
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FIG. 3: Fit to determine the branching fraction at E* = 1.548 GeV, corresponding to mx = 0. 
The data are the same as in Fig. 0(a) but with finer binning. The solid line is the total fitted 
spectrum, and the dashed line shows the background contribution to the fit. 



polynomial), range of the fit, and reasonable variations in event selection criteria. These 
systematic variations were added in quadrature to estimate the systematic uncertainty in the 
fitted yields. Systematic uncertainties were included in the branching fraction upper limits 
by scaling them upward by the same amount that Gaussian upper limits were changed when 
systematic and statistical uncertainties were added in quadrature. 

The resulting 90% confidence-level upper limits on the branching fraction are plotted in 
Fig. H] as the solid line. The dashed line shows the upper limit for statistical uncertain- 
ties alone. Statistical uncertainties are dominant for the full mass range in the plot, and 
systematic uncertainties are negligible for the low-mass points. The upper limits from this 
experiment are well below the previous experimental limits for charmonium radiative decay 
to invisible final states j22| . 

The present data yield a branching fraction for J /if} — > 7 + invisible of (1.5 ±2.4) x 10~ 6 
at m x = 0, and an upper limit of 4.3 x 10~ 6 at the 90% confidence level. The error quoted 
here (2.4 x 10~ 6 ) is from the statistical error in the fitted peak only. Substituting this upper 
limit in Eq. (1) gives cos 2 #/tan 2 /3 < 0.084. When combined with experimental limits for 
T radiative decay [ij}, following the procedure in Ref. Q, the present data suggest a very 
small value for cos 2 9. However a quantitative estimate cannot be made because only upper 
limits for the T branching fraction are available. Further measurements of T radiative decay 
to invisible states are warranted. 

In summary, the branching fraction for J /if) radiative decay to invisible particles has been 
measured as a function of particle mass. No statistically significant strength was observed 
for narrow states with mass less than 960 MeV/c 2 . The resulting upper limits place new 
constraints on some theoretical models that include a very light neutralino dark-matter 
candidate. 
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FIG. 4: The 90% confidence-level upper limits for J ftp — )■ 7X, where X is invisible to the detector. 
The dashed line shows the results for statistical uncertainties alone, and the solid line includes 
systematic and statistical uncertainties. 
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